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Abstract 

In this paper we study global Poincare inequalities on balls in a large class of 
sub-Ricmannian manifolds satisfying the generalized curvature dimension inequality 
introduced by F.Baudoin and N.Garofalo. As a corollary, we prove the uniqueness 
of solutions for the subelliptic heat equation. Our results apply in particular to CR 
Sasakian manifolds with Tanaka-Webster-Ricci curvature bounded from below and 
Carnot groups of step two. 
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1 Introduction 



In this paper, M is a C°° connected finite dimensional manifold endowed with a smooth 
measure p,. Let L be a second-order diffusion operator on M, which is symmetric, non- 
positive, locally subelliptic in the sense of [17], [25], [18], with LI = 0. 

By the subellipticy of L, there is an intrinsic distance d(x, y) associated to L on M 
([17], [25]). If we denote T(f) := T(f,f) (the Bakry-Emery's carre du champ of L, see [3]) 
by the quadratic differential form T(f,g) = \{L(fg) - fLg - gLf), f,g G C°°(M), the 
distance d(x, y) is defined via the notion of subunit curves with respect to the length of 
the gradient, yT(/). Throughout this paper, we assume that our sub-Riemannian metric 
space (M, d) is complete. 

This class of operators includes Hdrmander type operators of order k (or sum of squares 
of vector fields satisfying Hdrmander condition), Grushin's operator and sub-Laplacians 
on Sasakian CR-manifolds (see [25], [7]). For domains whose diameters are bounded in 
terms of the subellipticity constants, the Poincare inequality was proved in [24] for the 
Hdrmander type operator, and the general subelliptic Poincare inequality can be found in 
Kusuoka and Stroock's work [27] through their probabilistic methods. 

In general, the global Poincare inequality could fail if the domain is not bounded ([24], 
even in the Riemannian case). On the Riemannian manifold with Ricci curvature bounded 
from below by — K < 0, the global Poincare inequality is proved by Buser [12]: 

/ \f-f B{x , r) \ 2 dp<C ir 2 e c ^ r [ \Vffdp, V/eC°°(M),r>0. (1.1) 

JB(x,r) JB{x,r) 

This was done through geodesic arguments equipped with geometric tools such as the 
Bishop-Gromov comparison theorem. 

In our subelliptic framework, due to the lack of several properties such as ellipticity of 
L or smoothness of the distance function, we cannot use many of the tools for Rieman- 
nian manifolds. As a recent breakthrough, the generalized curvature dimension inequality 
CD(pi, p2, k, d) of Baudoin and Garofalo [7] was introduced to specify a certain subelliptic 
curvature condition ([7], [4], [5], [6], [8], [9], [10]). (Definition 2.1 in the present paper) 

For example, in [7] it was shown that if the lower bound of Tanaka-Webster-Ricci tensor 
on a CR Sasakian manifold of codimension 1 is given by p\ 6 K, the sub-Laplacian of M 
satisfies CD(p\, d/4, l,d) where the real dimension of the distribution of the CR manifold 
is d. 

If pi is non-negative (an analogue of the non-negative lower bound of Ricci tensor in 
the Riemannian manifold), some types of Poincare inequalities were discussed in [5], [10], 
including Buser's Poincare inequality with K = 0. 

Our purpose in the present paper is to prove the Buser's global Poincare inequality in 
the subelliptic framework with the condition CD(—K, p2, k, d), K > 0: 

Theorem 1.1 (Poincare inequality). IfM satisfies CD(—K,p2,K,d) with K > 0, for any 
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r > 0, x G M and f G C°°(M), 




where C p ±, C P 2 > depend on p2, k, d. 

Once the global subelliptic Poincare inequalities and the volume doubling property are as- 
sumed, one can obtain Sobolev inequalities through [36], and we can run Moser's iteration 
([31], [32]) to recover various properties such as the uniqueness of the solution for the heat 
equation through the mean value estimates of subsolution of L-Laplace equation. Note 
that two sided heat kernel bounds and parabolic Harnack inequality for P t f were already 
obtained through Li-Yau type inequalities and heat kernel methods in [7], [6]. However this 
Harnack inequality does not directly yield the uniqueness of positive solutions in contrast 
to Li-Yau's [29]. 

We consider the following heat equation associated with L. 



Adapting the ideas of [14], [30], we conclude our main corollaries, the uniqueness of the 
solutions for the subelliptic heat equation. 

Theorem 1.2 (Uniqueness of positive solutions). Assume that M satisfies the generalized 
curvature-dimension inequality CD(—K,p2,n,d), K > with respect to L. For any non- 
negative solution u G C(M x [0, oo)) of (1-3), u is uniquely determined by the initial data 
u(-,0) = f and u(x,t) = Ptf{x), where Pf is the heat semigroup generated by L. 

Theorem 1.3 (Uniqueness of L p solutions). For 1 < p < oo, L p solution u(x,t) of 

(1.3) defined on M x (0, oo) is uniquely determined by the initial data f G L P (M), i.e. 

LP 

u(-, t) G L P (M) for any t > and u — > f as t ->■ 0. 

If CD(—K, p2, k, d), K > is satisfied, we have an unique L 1 solution u(x,t) of (1.3) with 
the initial data u / G L 1 (M) as t — > 0. 

Theorem 1.2 improves Li-Yau/Harnack inequality for the subelliptic L which is proved 
for Ptf in [7], [6]. (Find notations in the following section. See Remark 2.5 for A e .) 

Corollary 1.4 (Li-Yau type inequality, [7]). Assume CD(pi, p2, k, d), p\ G R. Any non- 
negative solution u(x,t) G A e of (1.3) satisfies the following inequality: 



Corollary 1.5 (Harnack inequality, [6], [7]). Assume CD(—K,p2,K,d), K > 0. Then any 
non-negative solution u ^ of (1.3) satisfies for x, y G M, s < t G M+, 




(1.3) 




u(x, s) 
u(y,t) 



< 




exp f — (t - s) + —d(x, y) 2 + 



d(x,y) 2 
4(t - s) 
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2 Preliminaries 



Our subelliptic operator L is described as follows: L is a second-order diffusion operator 
with real C°° coefficients on M. There exists a neighborhood U of x G M and a constant 
C > 0, such that for any / G C§°(17), 

ll/ll? <C(|(/,L/>| + H/II1), (2.4) 

where ||/|| e = (/ |u(f)| 2 (l + \t\ 2 ) e d£) 1/2 is the Sobolev norm of order < e < f, and 
(■,■), || ■ || 2 are respectively the inner product and the norm of L 2 (M, y). Also L is sym- 
metric, non-positive and has zero order term, i.e.: 

/ fLgdfi= f gLfdfi, [ fLfdfi < 0, LI = 0, 

JM JM JM 

for every f,g G Cg°(M). 

The intrinsic sub-Riemannian metric associated with L is defined by the minimal length 
of subunit curve: 

d(x,y) = inf{T|3 Lipschitz 7 : [0, T] -»• M, 7 (0) = x,7(T) = y, 

l^/(7(*))l < Vr/(7(*)),V/ G C°°(M), almost every t G [0,T]}, 

where T(f,g) = ^(L(fg) — fLg — gLf), T(f) = T(f, /). We assume that the metric space 
(M, d) is complete. 

Note that this class of operators strictly includes Hdrmander's type operators 
— YlT^i X*X{ ( Xj's are the smooth vector fields satisfying Hormander condition of order 
k on M, and X* is the formal adjoint of Xi in L 2 (M, y) ). Following Strichartz [38], the 
completeness assumption of (M, d) yields that L is essentially self-adjoint on Cq°(M). So 
we can denote by L the unique self-adjoint extension (the Friedrichs extension) of L in 
L 2 (M, ju). Maximum principle([ll]) and Hormander 's hypoellipticity of L are well-known. 
See [17], [24], [7], [33] for more properties of L. 

We follow the steps in [7] to introduce the curvature assumption on our subelliptic frame- 
work. In addition to T, we assume that M is endowed with another smooth symmetric 
bilinear differential form, indicated with T z , satisfying for /, g G C°°(M) 

r z (fg,h) = fr z (g,h)+ g r z (f,h), 

and T z (f) = T z (f,f)>0. 

We make the following assumptions that will be in force throughout the paper: 
(H.l) There exists an increasing sequence hk G C^°(M) such that hk f~ 1 on M, and 

||r(/i fc )||oo + ||r z (/ife)||oo^0, as/c^oo. 
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(H.2) For any / G C°°(M) one has 

r(/,r z (/)) = r z (/,r(/)). 

(H.3) For every t > 0, P t l = 1 and for every / G Cg°(M) and T > 0, one has 

sup ||r(pj)|| 00 + ||r z (p t /)|| 00 <+oo, 

te[o,T] 

where p is the heat semigroup generated by L. 

(Details about the assumptions are discussed in [7]) The assumption (H.l) is implied 
by the completeness of the metric space. In the sub-Riemannian geometries covered by 
the present work, the assumption (H.2) means that the torsion of the sub-Riemannian 
connection is vertical (for instance, Sasakian condition of CR manifolds). Removing this 
assumption in certain cases is discussed in [9]. Assumption (H.3) is necessary to rigorously 
justify the Bakry-Emery type arguments. It is a consequence of the generalized curvature 
dimension inequality below in many examples (see [7]). 

In addition to V and T z , we denote the following second order differential bilinear forms: 
for tmy f,g€C°°(M), 

r 2 (/, g) = l [LT(f, g) - T(f, Lg) - T(g, Lf)] , (2.5) 

If (/,</) = \[LT Z U,9) -r Z (f,Lg) -T z (g,Lf)]. (2.6) 

As for T and T z , we denote r 2 (/) = T 2 (f, f), T z (f) = T z (f, f). 
The following curvature dimension condition was introduced in [7]. 

Definition 2.1 ([7], generalized curvature dimension inequality). We say that L satis- 
fies the generalized curvature dimension inequality CD(pi, p2, k, d) on M if there exist 
constants pi G R, p 2 > 0, k > 0, and < d < oo such that the inequality 

r 2 (/) + frf (/) > \{Lff + ( Pl - ^) r(/) + P2 v z (f) 

holds for every f G C°°(M) and every v > 0. 

The inequality CD(p±, p 2 , K,d) turns out to be equivalent to lower bounds on intrinsic 
curvature tensors in [7]. The following is an exemplary curvature condition implying 
CD(pi, p2, k, d) on CR manifold. 

Proposition 2.2. [7] Let (M, 9) be a complete CR Sasakian manifold with real dimension 
2n + 1. The Tanaka- Webster Ricci tensor satisfies the bound 

Ric x (v,w) > pi\v\ 2 ,\/x G M,Vv G H x , 

if and only if the curvature dimension inequality CD(pi, |, 1, d) holds with d = 2n and 
T z (f) = (Tf) 2 and the hypothesis (H.l), (H.2), (H.3) are satisfied. 
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With the curvature inequality condition assumed, various aspects on sub-Riemannian 
manifolds have been discovered in [7], [4], [5], [6], [8], [9], [10]. In particular, we have the fol- 
lowing essential properties - two-sided heat kernel bounds and volume doubling property 
of balls with exponential term. 

Proposition 2.3. [6] If we assume the curvature condition CD(—K,p2,n,d), K > on 
M, for any x, y G M, t > 0, r > 0, 

n(B(x, 2r)) < C dl eMC d 2Kr 2 )fi(B(x, r)), (2.9) 

where D = d ^1 + ^0 , Ci, C2, C3, C4, C^i, Crf2 are positive and determined by p2, k, d. 
Denote Q = log 2 C d \. (2.9) implies that for any A > 1, 

rut \ ^^ Rog 2 Al-i 
/x(£(x,Ar)) ^ri°g 2 Al f/ n, K ST rov^ 

V(B{x,r)) ^ (2J0) 

<C dl \QeM^K(\rf). 

This doubling property allows us to estimate /j,(B(x,y/t)) by fi(B(y,\/i)): 
fi(B{x, y/t)) < fJ,(B(y, Vi + d(x, y))) 

< »(B(y, Vt))2C dl (l + ^E^j ^ exp (^K{t + d(x, yf] 

So we modify the upper bound of heat kernel (2.8) with the volume of a single ball, i.e., 
for C5, Cq > depending on /?2, ft, d, 

Note that 1 + A < C(e)e eA for Ve > 0, A > is applied. 

Remark 2.4. As mentioned in [6], the square in the exponent of volume doubling prop- 
erty might not be optimal. For instance, in the Riemannian manifold with Ricci ten- 
sor bounded below by —K < 0, by the Bishop- Gromov comparison theorem we have 
V(x,\r) < V(x, r)\ n exp( yj (n — l)K(Xr)) where A > 1 and V(x,r) is the Riemannian 
measure of the ball B(x,r). 
This yields the difference of the exponent in (1.1) and (1.2). 
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Remark 2.5. ([4]) Notice that the positive solution u carries additional condition in 
the Li-Yau type inequality, Corollary 1.4- Due to the technical reason in the proof 
of Theorem 6.1 in [7], u needs to be contained in A t = {/ G C£°(M) : / — e > 
0, \/T(f — e), y/T z (f — e) G L 2 (M)}. Same restriction is required for log-Sobolev inequal- 
ity in [4j- 



3 Poincare inequality on the ball 

3.1 Lower bound of the Dirichlet heat kernel on the ball 

Throughout this section, L satisfies CD(—K,p2,K,d), K > on M. 
To adapt Kusuoka and Stroock's idea [27], the necessary ingredients will be two-sided 
heat kernel bound (2. 7), (2. 8) and doubling (2.10). 

Denote B = B(xo,r), sub-Riemannian ball centered at xo with radius r. On the ball B, 
the Dirichlet heat kernel pf ,D (x,y) will be defined by the transition probability 

p?> D (x,y)dp(y) = P[( > t,X(t) G d M (y)], 

where X(t) is the associated Markov process of the semigroup operator P t = e tL with 
X(0) = x, and the lifetime of X in B is C = mf{t > 0, X(t) B}. 

First, the lower bound of the Dirichlet heat kernel for close x, y can be obtained by the 
argument of Kusuoka and Stroock [27]: 

Lemma 3.1. For any k £ (0, 1), there exists C a = C(k, p2, n, d) £ (1, oo) such that for any 
xq G M, r > and a = c p^ -pjy G (0, 1), the Dirichlet heat kernel on B = B(xQ,r) 
has lower bound 

-<^^-(-c^) (, 12 , 

for all t G (0, (or) 2 ] and x,y G B(xo,kr) such that d(x,y) < ar. 
Here c, C > depend only on p2, k, d. 

Proof Let a = (C a (Kr 2 + l))" 5 G (0,1) with some C a > 1 which will be determined 
later. Note that K(ar) 2 < C^ 1 < 1 . 

Let d(x,y) < ar and t < (ar) 2 . The Dirichlet heat kernel can be written by the heat 
kernel of M and the lifetime of the process in the domain. That is, 

pf '°(x, y) = p t (x, y) - E x [ Pt ^(X(C), y),(<t], ( = mf{t > 0, X(t) B(x ,r)}. 

The lower bound (2.7) on the heat kernel pt(x,y) over the whole manifold yields 

C ie - 2C * ( Dd(x,y) 2 

Pt{x,y) >—— 7=77 exp 



■p(B(x,Vt)) "\ 2d t 
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If we use upper bound (2.8) on the heat kernel in the expectation, we have 

Pt-d x (0,y) < 



The balls can be replaced by concentric balls using the doubling property (2.10) as follows. 
ti(B(x,y/t)) ^ »(B(X(0,3r)) 



< 



»(B(X(0,Vt^)) -»(B(x(0,Vt^)) 

»(B(x,Vi)) < _p(B(y,2ar)) 



<C dl 



exp( — - — (2or) ). 



With d(X(Q,y) > r(l — k) and t — £ < t < (ar) 2 , the above controls imply that 
E x [p t - ( (X((),y),(<t} 



< — — exp 



H(B{x,y/t)) 



r 2 (l-kf 
10* 



+ 6C d2 Kr z ) - a exp ( -v-^ 



The last term came from W 



( 7(3r)(2ar) 



V v^c ) exp (" 



r 2 (l-fc) 2 ' 

' io(t-C) , 



holds if we choose C u > [j^^jj " , v < 



\Q/2 



< C u exp (— v-^i) , which 



(i-fc) 2 

20 



Combining these upper and lower estimates with t < (ar) 2 , d(x,y) < ar, 
pf' D (x,y) 

r< -2C 2 ( _ D d(x,y 

. we I 2d — r 

> — eV 



fi(B(x,y/t)) 

where C = C 3 C dl C v C^ l e iC <v +c * +2C * . 
Choose a small enough for 1 — C exp 
Then we conclude 



pf' D {x,y)> 



^ ( -^f-¥y) + SC d2 Kr 2 -%) 

Ax,y) 2 



> 



2 1 



exp — C 



n(B(x,Vt)) 

for all t < (ar) 2 , d(x,y) < ar, where c = 2~ l C\e~ 2C2 , C = ^ > depend on p 2 ,n,d. 
For instance, if we pick large C a = C(p 2 , k, d,k) > such as 



C a > max 



g + ln(2C 3 C dl C u C^e 3C ^+ c ^ 2C2 ) + 6C d2 (2d (I- k) 



v + 



(i-fc) 2 
10 



D 10 
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then our choice of 

° 2 = C a {Kr* + iy (3 - 13) 

satisfies the estimates above. 

□ 

Next step is the lower bound of Dirichlet heat kernel for any x, y in the smaller ball which 
is followed by the chain argument. Note that our lemma holds for any r > with the 
exponential square of radius, while the classic lemma holds only for < r < 1. 

Lemma 3.2. For any < k < 1 and < 5 < 1, there exists < c < 1, C > such that 
for any xo 6 M and r > 0, the Dirichlet heat kernel on the ball B = B(xo,r) has lower 
bound 

P ?> D (x,y)> CeM - CKr2) 



fi(B(x ,kr)) 
for all x, y G B(xq, kr) and 5r 2 < t < r 2 . 

Proof. Choosing a € (0,1) of (3.13) in the previous lemma, for all t < (ar) 2 , x,y G 
B(xo,kr), d(x,y) < ar, 

fj,(B(x,Vt)) t 

Now let x,y be any points in B(xQ,kr) and 5r 2 < t < r 2 . Set n = [16a -2 ], then 
16a -2 < n < YloT 2 . 

We choose {d}i=o,i,- ,2n C B(xQ,kr) such that 

Co = in = X , 6n = y, 
r ar 

d(£k,tk+i) < ~ < —r- 
n 4 

Let r = "k- Since Vt < X' if % G v 7 ^); then d{r] k ,Vk+i) < ar. 

By the previous lemma, 

J I'm, n, V 



Pr ij] k ,r] k+l ) >—— 7 =rrexp(-C 

K B {Vk,V T )) t 



^ cC d ^ exp{-C d2 K(ar) 2 ) , ^ d(r] k , r) k+1 ) 2 , 



And we see that 
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Observing 

pf' D (x,y)> [ ... / p?> D (x, m ) 

■ P?' D {m,m) • ■ ■ P B,D {r)2n-i,y)dr)i ■ ■ ■ drj2n-i 

we obtain 

1 / 4 \ 2n 

Doubling property (2.10) yields 

1 1 C dl 1 fc ( 3exp(-^Kr 2 ) 



n(B(x ,y/r)) fi(B(x ,r)) fi(B(x ,kr)) 
Also since cC^ 1 exp(-8C) < 1 and n < 17 a~ 2 = 17 C a {Kr 2 + 1) from (3.13), 

(cC£ ex P (- C (^ + 8 ) " Cd2K{arA \ 
> exp (-^- - (8C - In(cC^) + C d2 ) • 34C Q (Kr 2 + 1)) . 
This concludes our lemma 

fl>D c'exp(-C'i^) 

where 

c' = C" 1 ^ exp ("^ - (8C - ^(cC- 1 ) + Q 2 )(34C Q )) , 
C = ^f- + (8C - lnCcC^ 1 ) + C, 2 )(34C a ) 
are determined by p 2 , d> 5- D 



3.2 proof of Theorem 1.1 

In this section, we utilize Dirichlet, Neumann heat semigroup which can be found in 
[40], [35], [27], [22], then we will follow the arguments in [27] to prove Poincare inequality 
(1-2). \ ' 

Let B = B(xq, r). Define a subspace D°° C C co {B) as a collection of functions / satisfying 
-f B gLfdfi = J B T(g,f)dn for V# G C°°(B). Note that Cg°(B) C D°° C C°°{B). 
The Dirichlet form £(f,g) = f B T(f,g)dfi on D°° is closable in L 2 (B), and by closing 
it we gain a Dirichlet form and associated Markov heat semigroup P t ' with Neumann 
boundary condition. 
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B TV 

If we denote p t ' by the Neumann heat kernel over B, it will be a smooth kernel of the 
Neumann heat semigroup and its associated transition probability function. Naturally, 

since C^°(B) C D°°, the Neumann heat kernel dominates the Dirichlet heat kernel, i.e., 

B.N ^ B,D 
Pt >P t ■ 

Proof of Theorem 1.1 . We will prove the inequality with B(xQ,r/2) on the left hand 
side. Then by the Whitney type covering lemma (section 5 in [24]), we can match the 
balls on the both sides. The Whitney decomposition only requires a doubling property in 
the domain of argument. In B(xq, lOr), the doubling property holds with fixed constant 
Cdi exp(C(i2-f^(10r) 2 ), which will be multiplied at the end following the argument. 
From the previous lemma, for x,y G B(xo,r /2), 



P% M ' N ( X ,y)> 



ce -CKr* 



fi(B(x ,r/2)Y 
For any / G C°°(B) and x G B(x ,r/2), 

P* M > N U-P* {X0 > r) ' N f{x))\x) 



-• / (f(y) - P* iX0 ' r) ' N 'f{x))*Mv) 

, r / 2 )) JB(x ,r/2) 



ce 

n(B(x , 



'B(x ,r/2) 

-CKr 2 

( < 

> 



ce~ u f 

H{B(xo,r/2)) J B (x ,r/2) 

On the other hand, 

/ P^' r) *{f-P^' r) *f{x)ftz)M*) < I (f 2 - P r B 2 M > N f(xf)dn(x) 

JB(x ,r) JB{x ,r) 

JB(x ,r) al 



f I 

Jo Je 



-2Pf M ' N f(x)LP t B{xo ' r) > N f(x)dv(x)dt 

B(xo,r) 



= r I 2T(P t B ^ N f(x))d^x)dt 

JO JB(x ,r) 

< 2r 2 I r(/)d/i, 
JB(x ,r) 

where the last inequality comes from iT{P t f) < 0. And we obtain our desired conclusion 
/ (/(*) - fB(x ,r/2)) 2 dKx) < C' pl r 2 e C ^ Kr2 [ r(/)d/i, 

JB(x ,r/2) JB(x ,r) 

with C' pl = 2/c, C' p2 = C. □ 
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3.3 Sobolev inequality and L p mean value estimate 

As a consequence of Theorem 1.1, this section is dedicated to the Sobolev inequality 
and L p mean value inequality for subharmonic functions, which will be essential to prove 
Theorem 1.2 and 1.3. Throughout this paper, harmonic (resp. subharmonic) functions are 
/ G Dom(L) satisfying Lf = (resp.L/ > 0). 

Assume that L satisfies CD(—K,p2,K,d), K > on M. We have Poincare inequality 
(1.2) and exponential doubling property (2.9). 

With these two ingredients, one can derive local Sobolev inequality in [36], [37]. This is 
a classic path to Moser's iteration for Harnack's inequality. See theorem 2.2 in [36] and 
section. 10 in [37] (also the last section of [42]). 

Note that in [19], one can find that the weak L 1 Poincare inequality with the dou- 
bling property derives the isoperimetric inequality (and Sobolev inequalities) in Carnot- 
Caratheodory spaces. 

Proposition 3.3 ([36], [37], Sobolev inequality on balls). If the Poincare inequality (1.2) 
and the volume doubling condition (2.9) are satisfied for any r > 0, then for any x G M, 
< r, f G Co°(B(x,r)), denoting B = B(x,r), 

(jm L l/|ft<i ") V £ CirfC ' K " (m L im + r " l/|2)<i ") ' • (3 ' I4) 

where Q = log 2 Cdi in (2.9), C\,C e > depend only on p2,n,d. 
Note that this Sobolev inequality can also be obtained by following steps of [37]. 

With CD(—K,p2,K,d) assumed, by the upper bound of the heat kernel (2.11) 

pf >°(x, y) < p t (x, y) < ^P (<hK{t + d(x, yf) - 

where B = B(xo,r). Since < t < r 2 and d(x,y) < 2r for x,y G B, by (2.10) 

/ \ Q 

p(B(x ,r)) < p(B(x, 2r)) < C ( -^=) exp(c^r 2 )^(5(x, Vt)). 
Therefore, the Dirichlet heat kernel will be bounded from above by 

p?> D ( X ,y) < ( L^ e ^r- < C> Q^-Q/2 c'Kr 2 

Pt K ,y > ~ fi(B(x,Vt)) ~ p(B(x ,r)) 
Proposition 10.1 in [37] (also [42]) states that 

\\P t B ' D h^oo <C r Q/2 , V0<t<t 

< Cl /Q (C\\y/f(jj\\ 2 + 1/2 ||/|| 2 ) , V/ G C^(B). 
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Taking C = ^ B % 0<r) f Q e dKr2 and t = r 2 with A + B < (2A 2 + 2B 2 ) 1 / 2 , the Sobolev 
inequality (3.14) is proved. 

Once the local Sobolev embedding is acquired, our goal of this section, L p mean value 
estimate, can be obtained through the Moser's iteration. One can find the arguments for 
the Riemannian case in [35]. 

In [37], [36], the author obtained parabolic L p mean value estimate. But in our context, 
LP mean value estimate for subsolution of L-Laplace equation will be enough. 

Lemma 3.4 (One step of Moser's iteration). We assume that M satisfies 
CD(—K,p2,K,d), K > 0. For any subharmonic function u(x) > 0, i.e. Lu{x) > ; 
andO < Ri < R 2 < R, p> 2, 



[ u p °d»<C 2 e 2C ^—^—V^(\ u p d») , (3.15) 



where = 1 + %, B(-) = B(x , ■), V = n(B(R)). 

Remark 3.5. For any < R\ < R 2 < oo, there exists a Lipschitz continuous cut- 
off function t[ > satisfying i[\b(r 1 ) = 1, suppV' C B{R 2 ), y/T(ip) < R ,^L Rl almost 
everywhere for some C > which is independent to R\,R 2 . See theorem 1.5 in [20], 
lemma 3.6 in [13] and [39]. 

Proof. Denote B = B{R). Let ijj > be a cut-off function satisfying V'ls^i) = 1> supp ip C 

b(r 2 ), yfxvo < 

ij 2 ^i? 1 almost everywhere. Choose a test function <p — ip 2 u p 1 > 0, and 
the subharmonicity condition of u implies 

< (Lu, <p)= [ -T(u, ^ 2 u p - l )dp = [ (-(p - l)t[ 2 u p ~ 2 T(u) - 2^ p_1 r(n, ^)) dfi. 

JB JB 
Hence, by Cauchy-Schwarz inequality 

P~ 



- / ^ 2 u p - 2 v{u)dp< I Vu p_ Vr(u)r(V0<^ 

JB JB 

<(J t[ 2 u p ~ 2 r(u)dp^j 1 (J u p r(ij)df?j 1 

So, we have 



(p-1) 5 



f ^ 2 u p - 2 T(u)dfi < [ u p T(i[)dp. (3.16) 

JB JB 



4 

On the other hand, if we apply Holder inequality and the local Sobolev inequality (3.14) 
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li{B) 
< C 2 R 2 e 2C * KR2 



< CiR 2 e 2C * KR2 



Q-2 _2 



on ifju p / 2 , it will give 



_2_ 

_4 / 1 /• \ Q 

Q 
oo 



'supp(^) 

Using (3.16), the gradient term can be written by 

J T(4m p/2 )dv < J 2{u v T{^) + ^ 2 u p - 2 T{u)^ dfi 



2p 2 

< 2 + 



(P " I) 5 



2 



Therefore, given suppV> C -B(i? 2 ), V'ls^i) = 1, < ^ < 1 and ||r(-0) ||oo < { r 2 - Ri )^ we 
obtain 



— --— [ u pe dfi 
M- 8 ) JB{Ri) 

2p 2 \ C 2 R 2 



< C 2 e 2C * KR2 ( ( 2 + 



(p-i) 2 j (R 2 -R 

where 6 = 1 + ^. Note that we can assume that C > 1 without loss of generality. 

The desired inequality (3.15) is proved with C2 = \\C 2 C 2 . □ 

Now by iterating the above lemma, we prove L p mean value estimate. 

Theorem 3.6 (LP mean value inequality, p > 2). For any < S < 1, any p > 2, and any 

non-negative subsolution u of Lu = in a ball B(R) of volume V , 

supK} < C 3 e QCsKR2 (l - 5)- Q (v- 1 [ u p dfi) . (3.17) 

SB \ JB J 
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Proof. For i = 0, 1, 2, ■ ■ ■ , set pi = p6 l where 6 = 1 + ^. 
And let R = R, Ri - R i+ i = %+\f , i.e. 



3=1 

By Lemma 3.4, 



JB{R l+1 ) I 1 _ d ) \JB{Ri) J 



This yields 



Simple computation shows that 



oo 



l i Q j Q(Q + 2) 

3=1 3 =1 



1 



lim 

i—>oo 



r \ p«+i 

/ u Pi+1 d[i = sup{it}, 

JB{R l+1 ) I SB 



Q. Q(Q+2) 

Conclusively, where C3 = C 2 2 2 , we have 



sup{u} < (c 3 e QC * KR \l - d^V- 1 )* ( [ u p dA " . 

SB ^ ' \JB(R) J 



□ 



Corollary 3.7 {L p mean value inequality, < p < 2). L p mean value inequality (3.17) 
also holds for any < p < 2 with the constant C3 replaced by some C4 = C(Q,p). In 
particular, for p = 1 

sup{u} < C m e c ™ KR \l - S)~ Q f-^r / 
sb Je 

where C m , c m > depend only on P2, k, d. 



udjji ) , 

B 
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Proof. Let \ < a < 1 and p = a + (1 - a)/4. 

By Theorem 3.6 (we can pick Rq = pR in the proof) and J B u 2 dp < sup B {u 2 ~ p } J B u p dp, 
sup{n} < Ce cKR \l - <t)-% (v 1 I u 2 dp) * 

oB \ JpB J 

< (cV~ 1 2{ [ u p dp)h cKR2 ) (l-o-)- 2 supju 1 -!}. 

V JB J pB 

Now fix 5 G (§, 1) and set a = S, cr i+1 = a { + (1 - (Jj)/4 = cr, + (§)*(1 - <5)/4. Then 

sup{u} < A(^) Qi / 2 (l - sup {n}) 1 "!, 

o-jS 3 a i+1 B 

where A = (c^ 1 f B u p dp)h cKR2 ^j . 
Finally, the same iteration of Theorem 3.6 yields 

4 2Q 2 _Q 

sup{n} < (-) p 1 Ap (1 — 6) p 
SB 3 



42Q2 2c 2 Q / f \ p 

(-)^c?e^ KR (1 - 5)~T (V- 1 J u p dp\ . 



□ 



4 Uniqueness of the positive solution 

4.1 Minimality of the heat semigroup for positive solutions 

To prove Theorem 1.2, we reduce the question to the zero initial data. Following Lemma 
4.1 enables the reduction. This section is based on the idea of [14]. 

Lemma 4.1 (minimality of the heat semigroup). Let u 6 C(M x (0, T)) be a non-negative 
super solution of the heat equation (1.3) with initial data f £ Lf oc (M.), f > 0. 

i? 

Then P t f(x) = J M Pt(x , y) f (y)dp(y) is a smooth solution of (1.3) satisfying P t f — ^» / as 
t^O and u(-,t) > P t f. 

Proof. For any SI <s M, we denote P t ,D the Dirichlet heat semigroup associated with L 
on SI. Using the maximum principle for P t ' f — u(-,t), we have 

u{x,t)> P^ D f{x), VxGfi 

Denote fk = /ln fc £ L 2 (M) for the exhaustion {Sl^}. As shown above, u(-,t) > p^ k > D f > 

P^ k,D /j for all i Since P^ k ' D f i P t f i as — >■ oo, we have u(-,i) > P^/i almost 

everywhere for all i. Therefore, 

u{-,t) > Ptf almost everywhere. 
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To prove that the smooth P t f solves the heat equation, first we see that Ptf G L ; 1 oc (M) 
from the above estimate. Denote 

u k = P t (min(/, k)ln k ), 

then Uk is a smooth solution of the subelliptic heat equation and Uk /• P t f as k — > oo at 
any (x,t) G M x (0,T). 

For any if G C °°(M x (0,T)), since P t / G L[ oc (M), 

|(%> + L<p)u k \ < (sup |%> + L^\)l snppip P t f G L^M), V* G N. 
This allows us to take the limit of the integrand on the left hand side of 

/ / (dfip + L(p)ukd[idt = / ip (L — dt) Ukdpidt = 0. 
Jo Jm Jo Jm 

Therefore Ptf is a distributional solution of the subelliptic heat equation, and also it is 
smooth by the smooth convergence of u& to P t f and the hypoellipticity of L — dt- 
Once the smoothness of Ptf and u > Ptf are proved, the initial condition is straightfor- 
ward as follows : On any Q <e M, 

Pt(fln)<Ptf <«(-,*)■ 
When t -> 0, u ^% f and P t (/l n ) ^> /l n . Hence P t f f. □ 

4.2 proof of Theorem 1.2 

From the minimality Lemma 4.1, for any non-negative continuous solution u of (1.3), 

w(x, t) = u(x, t) — Ptu(x, 0) 

is a non- negative solution of (1.3) with zero initial data. Thus we can reduce the uniqueness 
of the positive solution to the zero initial data case. 

Let w(x, t) be any non-negative solution of the heat equation (1.3) with initial data / = 0. 
Define v(x,t) = fQw(x,s)ds. Our goal is to show v = 0, and so is w. 

Remark 4.2. v(x,t) = J * w(x, s)ds is a non-negative solution of the heat equation (1.3) 
with zero initial data, and subharmonic in x, i.e. Lv(-,t) = f Q l Lw(-, s)ds = w(-,t) > 0. 

The following growth estimate condition is originally suggested by Tikhonov for the 
uniqueness of the solution for the heat equation. 

Proposition 4.3 (Growth estimate of the solution, Tikhonov's condition). For any e > 

and < t < e, if v G C(M x (0, e)) is a non-negative solution of the subelliptic heat 
equation (1.3) satisfying Lv(-,t) > 0, then 

v(x,t) < Ciexp(C 2 d 2 (p,x)), 

where C\ = Ci(e) > 0, C2 = C^e) > 0, and d(p, •) is the distance from a fixed p G M. 
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Proof. Let B = B(x,d(p,x) + 1). Fix T > 0. From the minimality Lemma 4.1, 

+ T) > P T v(-,t) = / pr(p,y)v(y,t)dfi(y) > / p T (p,y)v(y,t)dfj,(y). 

JM JB 

From the curvature condition CD(—K, p 2 , k, d), the lower bound of heat kernel (2.7) is 

Pr(p,y) > C 3 exp(-C 4 d 2 {p,y)), 

where C 3 = C 3 (p, T, K, p 2 , k, d) > 0, C 4 = C 4 (T, K, p 2 , k, d) > 0. 
By the triangle inequality d(p, y) < 2d(p, x) + 1 for y 6 B, 

[ v(y,t)dn(y) < C 5 exp(C 6 d 2 (p,x)) v(p,t + T). 

JB 

By L 1 mean value estimate of Corollary 3.7 for the subharmonic function v(-,t) 

v(x,t) < C 7 exp(C 8 d(p,x) 2 ) / v(y,t)dfj,(y), 

JB 

where C7, C$ > depend on K, p 2 , k, d. Therefore, we obtain 

t) < C 9 eMCiod 2 (p, x))v(p, t + T), 

where the constants depend on p, T, K, p 2 , k, d. As t varies from to e, v(p, t + T) remains 
uniformly bounded in t. So we have the desired conclusion. □ 

Together with the previous proposition, the proof of Theorem 1.2 is finished by the 
following proposition. 

Proposition 4.4. If v(x,t) is a solution of (1.3) with initial f(x) = satisfying 

\v(x,t)\ < CiexpC 2 d 2 (p,x) 
for some positive C\, C 2 , then v = 0. 

Existence of Lipschitz cut-off function and integration by part allow us to follow exactly 
the same proof of corollary 11.10 in [21]. 

5 Uniqueness of L p solution 
5.1 proof of Theorem 1.3, p > 1 

For p = 00, the uniqueness of the L°° solution, or equivalently the stochastic completeness 
of M, can be found in [7]. If p G (1, 00), without any curvature assumption, the uniqueness 
follows immediately by adapting the idea of [30]. 



18 



Theorem 5.1. Let v(x,t) is a non-negative function defined on M x (0, T) with 

d \ 

¥t -L)v(x,t)<0 

v as t -> 
v{-,t)eL p (M) V*€(0,T), 

i/ten i) = on M x (0, T) . 

In particular, any L p solution of the heat equation is uniquely determined by its initial 
data in L P {M). 

Proof. Fix xq G M an arbitrary base point. 

From remark 3.5, we choose if)(x) G Co(5(xo, 2i?)) a cut-off function satisfying i/)\b(x ,R) = 
1, < V < 1, Hv^Wlloo < § for some C > 0. 

Since u is a subsolution with the zero initial data, for any r G (0, T), 



/ (\p 2 (x)v p ~ 1 (x,t)Lv(x,t)dfi(x)dt> ( ( %l) 2 (x)v p ~ l ^dn{x)dt 
Jo Jm Jo Jm ot 

= ~ r W + ( I ^(xVdfiix)) dt = - [ i; 2 (x)v p (x,T)dn(x). 

P JO at \JM J P JM 

On the other hand, integrating by parts yields 
i) 2 {x)v p ~ l (x, t)Lv(x, t)dfi(x)dt 



2^v p - 1 T(^,v)dndt- / V 2 (P- l)v p ~ 2 T(v)dfidt. 
o Jm Jo Jm 



On the right hand side, observing 



we obtain the following estimate. 

^ 2 (x)v p ~ 1 (x, t)Lv(x, t)dn(x)dt 

i 

< r f -^—Ti^d^dt - r i ^-^ 2 v p - 2 v{ v )d^dt 

Jo JMP-1 Jo Jm 2 

= r [ -^—T^)v p d^dt - 2{p 7 x) r / ^ 2 T(v p / 2 )dfidt. 

io JmP-1 P Jo Jm 

Combining with the previous conclusion and the assumption | y/T(ip)\ < ^ , 

f ^ 2 (x)v p (x,r)dfi{x) + 2(P ~ ^ / T / tp 2 T(v p / 2 )dfidt < . 2pC ,l n0 f T I v p dndt. 
Jm P Jo Jm (p ~ 1)# 2 Jo 
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As R — > oo, since T(v p ^ 2 ) > 0, we have 

f v p (x,T)dfi(x) = Vt€(0,T). 

JM 

Thus, v = 0. □ 
5.2 Hamilton's inequality 

Before we move on to L 1 solutions, we will prove the gradient estimate of the logarithm of 
the heat kernel. We will apply subelliptic version of Hamilton's inequality which was orig- 
inally proved for closed Riemannian manifolds in [23], then for non-compact Riemannian 
manifolds in [26]. 

Proposition 5.2 (Hamilton's inequality). Assume thatM. satisfies the curvature condition 
CD(—K,p2,K,d). If a positive solution u £ A t to the subelliptic heat equation satisfies 
u < M on M x (0, T) for some M > and < T < oo, one has 

t r ( ln„(x,<))<(l + ^ + 2«)m(-^ y ) (5.18) 

for all (x,t) € M X (0,T). 

Proof. By Theorem 1.2, it suffices to show that the estimate holds for u = Ptf G A e > 0. 
(See Remark 2.5 for A e .) We apply the reverse log-Sobolev inequality in [4], i.e., 

tPtfT (In P t f) + p 2 t 2 PtfT z (In P t f) < (l + ^+2Kt) (P t (flnf)-(P t f)lnP t f). 

Then our desired inequality is instantly deduced by P t (f\nf) < (P t /)lnM and 
P2 t 2 P t fT z > 0. □ 

Lemma 5.3. // M satisfies CD(—K,p2,K,d), there exists Ch = Ch(p2,K,d) > 0, t > 
0,x,y G M, 

T x (ln Pt (x, y))<°±( 1 + *! + K t) (K(t + d(x, y) 2 ) + ^ vY 



t V P2 ) V * 

Proof. Let i > and y G M. Let u(x,s) := pt +s (x,y), then u is a smooth, positive 

solution to the heat equation. From the heat kernel upper bound (2.11), for to = qc 6 k -> 
t 

2' 



o < t < t , < s < |, Vx e M, 



u(x, S ) < Cs r exp (c 6 tf(£ + a + d(x, y) 2 ) - §^ 



»(B(y,Jl7s)) V 2 6(- 2+ s) 



< < = M. 

p(B(y,Jt 2 +s)) p(B(y, x n)) 
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Moreover u(x, s) < M for all s > 0, since ||Pt||oo-K» < 1 for any t > 0. 
By the Hamilton's inequality (5.18), the heat kernel lower bound (2.7) for u(x, s) with 
s = | and the doubling property (2.10), 

L„ . i xx / 2k \, / M 
-T(\nu(x,-)) <[l + — +Kt)\n' 



P2 J \U(X,%) 

<[i + ^ + Kt)^( K (t + d (x,yf) + ^l 



P2 J 2 

where C h = 2 In {C^ x C dl 2 Q l 2 ) max(^, ^ + C 2 ). □ 

If we combine the previous lemma with (2.11), we obtain the following simpler statement 
for small t, which will be useful in the next section. 

Lemma 5.4. Assume CD(—K,p2,n,d). For any R > 0, ft > and xq G M, there exists 
C > 0, t > such that for d(x, y) > R/A, <t <t Q , 



^T y ( Pt (x,y)) 



Ce~P R2 

< 



fi(B(x,Vt))' 



5.3 proof of Theorem 1.3, p = 1 

Prior to the uniqueness of L 1 solution for the heat equation, we prove the uniqueness of 
L 1 harmonic function. Basic idea of the proof comes from [30]. 

Remark 5.5. We assume the fixed curvature bound p\ = —K instead of the negative 
quadratic lower bound of Ricci curvature of [30]. 

Theorem 5.6. // M satisfies CD(—K,p2,n,d), then any L 1 non-negative subharmonic 
function on M must be identically constant. 

In particular, any L 1 harmonic function on M must be identically constant. 

Proof. Let g € L X (M) be a non-negative function satisfying Lg > 0, i.e. subharmonic. For 
any t > 0, 



LP t g(x)= (L x p t (x,y))g(y)dn(y) 

JM 

= / (-^:Pt{x,y))g(y)diJ,(y) = / (L y p t (x,y))g(y)dp{y). 

JM ol JM 

We claim the following integration by parts. 

/ (L y pt(x,y))g(y)dn(y) = / p t {x,y)Lg(y)d/j,{y). (5.19) 

JM JM 
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M 



< 



To justify the claim, we observe the following 

My) [g{y)L y Pt{x, y) - p t {x, y)Lg(y)] dfi(y) 

■ [T(4>Rg,Pt(x, •)) - r(if) R pt(x, •),£?)] dn 

- [gF(4>R,Pt{x, •)) - Pt(x, -) r (^H, 9)] dy. 

C (gy/T(pt{x, ■)) +Pt(x, ■) V^(i)) d/i. (5-20) 

lB(x ,R+l)\B(x ,R) v 7 

where V-R > is a Lipschitz continuous cut-off function satisfying V'i?l_B(x ,_R) = 1) 
suppV'i? C B(xq,R + 1) and yT'(V'-R) < (7 almost everywhere for some C > which 
is independent to R > 0.(See Remark 3.5.) 

It suffices to show that both integrals on the right-hand side vanish as R — > 00. We can 
consider R large enough so that x € B(xq,R/4). 

Let be a cut-off function for an annulus satisfying ^>\b(x () ,r+i)\b(x () ,r) = 1> 
¥ ? ls(x ,il-i)u(M\B(xo,il+2)) = and \/ r (^) ^ ^ alm ost everywhere. By the subharmonicity 
of g, 

0< (p 2 gLgdfi = -2 tpgT(<p, g)dfi - (p 2 T(g)d^ 
jm jm jm 

< [ \-h2g^¥g)-vVW)) 2 + ^(v)g 2 -\v 2 ng)\ ^ 

Jm L z z 
<2 / r^d/i-i / <p 2 T{g)d^ 

JM z JM 

Therefore, applying L 1 mean value estimate Corollary 3.7 to g, 

[ T(g)dfi < 4 f T(ip)g 2 dfi < AC 2 f g 2 d^ 

J B(x ,R+l)\B(x ,R) JM JB(x ,R+2) 



<4C 2 \\g\\ Ll sup g(y)<C'e c ' KR2 1 

B(x ,R+2) n(B(x ,2R + 4)) 



By Schwarz inequality, 

/ Vmdfi < ( [ 

JB(xn,R+l)\B(xn,R) \JE 



B(x ,R+l)\B(x ,R) 



r(^)d/i (/i(B(x ,212 + 4)))3 



Li- 



In addition to this estimate, to bound the second integration in (5.20) we consider the 
upper bound of heat kernel (2.11): 



Pt(x,y) < 



C 5 



H(B(x,y/t)) 



e W (C 6 K(t + d(x,y) 2 ) 



d{x,yf 
6t 



22 



Combining the above two inequalities, we estimate the second term of (5.20) for small 
< t < T = T(K,p 2 ,K,d). 

Pt(x, -)\/T(g)dfi 

B(x ,R+l)\B(x ,R) 



< sup Pt(x,y)\ / V T (a) d V 

\yeB(xo,R+l)\B(x ,R) J J B(x ,R+l)\B(x ,R) 

C ( R2 \ II II *-"*> n 



with d(x , x) < f, f < d(x, y) < AR and a > 0. 



For the first term of (5.20), L 1 mean value estimate for g yields 

/ gV^iMx^yjdfi 

J B(xt),R+l)\B(x(),R) 



< sup g) ^T(p t (x,-))dfi 

\B(xo,R+l)\B(xo,R) J Jb(xo,R+1)\B(x ,R) 

' Ce cKR2 
Ii(B(xq,2R + 2)) u:ju1j J Jb(xo,R+i)\B(x ,R) 



I O c J v J v \ f 

^ l —TFT, nn I oXN llffllz. 1 / y/F(pt(x,-))dfi. 

J JB(x ,l 



If we apply Lemma 5.4 for (3 > cK, 

< / p(B(x ,R + l)) \ Ce-?' R2 < Ce~^ 2 g-^ 

where < t < T = T(K, p 2 , k, d) small enough and /?' > 0. 

Therefore, as R — > oo, the integration of (5.20) vanishes as we desired, and we proved 
our claim (5.19). 

Now since the integration by part (5.19) holds for small t, we have 

j t P t g = LP t g = P t (Lg) > 0. 

And by the semigroup property, Ptg(x) > g(x) for all t > 0, x 6 M. 

On the other hand, by the stochastic completeness ([7]) of Pt, \\Ptg\lL 1 = WqWl 1 - 
Therefore Ptg = g, i.e. g is harmonic. 

For any constant 7 > 0, (5—7)+ = max(0, <?— 7) < g is also a non-negative L 1 subharmonic 
function. And by the same argument, it is harmonic. min(<7, 7) = g — (g — 7)+ is also 
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non-negative L 1 harmonic function. Observe that min(<7, 7) G C°°(M) for any 7 > by 
the hypoellipticity of L. This is not possible unless g is constant. 

Finally, any harmonic function u G L 1 (M) is identically constant since \u\ is non-negative 
L 1 subharmonic function which must be constant by the above. □ 

With the uniqueness of L 1 harmonic function, we are ready to prove L 1 uniqueness of the 
solution for the subelliptic heat equation. 

Theorem 5.7. Let M satisfy CD(—K, p2, k, d). Let v : M x [0, 00) — > E be a non-negative 
function satisfying 

- v(x, t) > 0, \\v(-,t)\\ L i {M) < 00, Vt > 0, 

l|w(->*)IUi(M) o> 

f/ien v(x,t) = on M x (0, 00) . 
Proof. For any e > 0, denote 

t) = max (0, t + e) — Pt(u(-, e))) . 

Then it follows that ^ e > 0, lim t ^ ipe(x, t) = 0, (L - J^) ip e > 0. 
Fix T > 0. Define 

f(x) = I ip e (x,t)dt, 
Jo 

which satisfies Lf{x) = ?/> e (x,T) — ij) € (x, 0) > 0. 

The assumption t>(-,t) G L 1 (M) yields J J M t + e)|d/x(x)dt < 00. Together with 

Jo Im P t v i x ^ € ) d ^ dt - T '/m^^M 1 ) < 00 ) we obtain ||/|| L i (M ) < 00. 
Now / is non-negative L 1 subharmonic function, so that we can apply Theorem 5.6 to 
/ and conclude / is identically constant. This implies = Lf(-) = ip e (-,T) for arbitrary 
T > 0. Hence for any t > 0, 

v(x,t + e) < P t (v(-,e))(x) 

<||p t (x,-)||oo||<,e)||Li<M||<,e)|| Ll ^0, 

where the uniform bound for ^Woo is found in Lemma 5.3. 

Therefore non-negative v(x,t) must be zero for all (x,t) G M x (0, 00). □ 

Proof of Theorem 1.3, p = 1. For any L 1 solution u of (L — t|) u = with the initial 
condition u / G L 1 (M) as i -> 0, 

v(x,i) := |u(ar,t) - P t /(x)| 

will be a non-negative L 1 subsolution of the heat equation with v as t — > 0. 
By the previous theorem, d e on M x (0, 00). Therefore, u is uniquely determined to 
be P t f. 

□ 
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Remark 5.8. In [2],[1], the measure contractive definition of Ricci tensor bound and vol- 
ume comparison theorem (which is not yet established in our framework) are introduced in 
three dimensional sub-Riemannian spaces. This measure contraction property is extended 
to higher dimensions in [28]. 

One can find the uniqueness theorem of the positive solution in symmetric local Dirichlet 
spaces, provided a local parabolic Harnack inequality in [16]. 

Acknowledgements. The author is grateful to his advisor, Fabrice Baudoin, for sug- 
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